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Abstract 

Lattice results, kinematical constraints and QCD dispersion relations are combined for the first time to 



' derive model-independent bounds for QCD form factors and corresponding rates. To take into account the 
' error bars on the lattice results we develop a general formalism which ascribes well-defined statistical properties 
to these bounds. We concentrate on B° — ► 7r + f~i^ decays because of the relative simplicity of the analytical 
■ behavior of the relevant polarization functions and because of their immediate phcnomcnological relevance. To 
determine the range of applicability of the dispersive matching required to obtain the bounds, we have evaluated 
the leading perturbative and non-perturbative QCD corrections to the relevant polarization functions. Despite 
the lattice results' large error bars and limited kinematical range, the bounds we obtain appear to favor certain 
parametrizations for the form factors over others. The bounds also enable the determination of \V u b\ from the 
total experimental rate with a theoretical error ranging from 27% to 37%, depending on the assumptions made. 
The techniques developed here are, in fact, quite general and are not limited to use with lattice results nor to 
p> ■ semileptonic B° — > ir + £~i>i decays. 
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1 Introduction 



The CLEO Collaboration has very recently reported preliminary measurements of the total rates for semileptonic 
B° — > tt + £~Di and B° — > p + £~Pg decays m. ^] The study of these rates is important, for it will eventually lead 
to accurate determinations of the magnitude of the poorly known Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element V u b- f] Such determinations, however, require precise calculations of the relevant form factors which 
are non-perturbative QCD functions of q 2 , where q is the four-momentum transferred to the leptons. These 
calculations are difficult because they involve understanding the underlying QCD dynamics over a large range 
of momentum transfers from q 2 = q^ax — ( m B ~ fn-n( P )) 2 = 26.4 GeV 2 (20.3 GeV 2 ), where the final state hadron 
is at rest in the frame of the B meson, to q 2 ~ 0, where it recoils very strongly. 

There exist several quark model determinations of these form factors. The authors of Ref. || (ISGW) 
determine these form factors around q 2 nax in a non-relativistic model and extrapolate them in a somewhat ad 
hoc manner to the fully relativistic region around q 2 = 0. The WSB model Q is a first attempt at a relativistic, 
quark-model treatment of these decays, but only determines the form factors at q 2 = and extends them to 
arbitrary q 2 by assuming pole behavior. The authors of Ref. || follow the WSB method but use pole or dipole 
behavior according to the power-counting rules of QCD. Since these rules do not modify the WSB results for 
B° — > TT + £~De decays, we shall not consider the results of Ref. Q any further. The authors of Ref. || also obtain 
the relevant form factors at q 2 = 0, in a light-front quark model. In Ref. j?J these form factors are calculated in 
the ISGW2 model which is a version of the ISGW model that is more consistent with Heavy Quark Symmetry 
(HQS) constraints, more consistent with relativity and better behaved at large recoils. None of these quark 
model calculations, however, determines the full dynamics of the decay process. Only very recently have light- 
front quark models begun to be used to determine the g 2 -dependence of B — > iriP^ decay form factors over a 
large range of q 2 jg, |s) . 

The form factors relevant to B — > tt£Di decays have also been calculated using QCD sumrules. Most of these 
calculations determine the form factors at one or two values of q 2 and assume a g 2 -dependence to obtain the 
form factors at arbitrary q 2 fl0||-||i3||. In Ref. jlj|, though, this g 2 -depedence is calculated to leading order in 
and in Refs. |l4| ] and |l5| over a relatively large rang e of q 2 from to about 15 - 20 GeV 2 . The authors 
of Ref. jl5| even suggest, in Ref. ]l6| , that their result for the form factor / + (<? 2 ), which dominates the rate for 
B — > ir£i?i decays, could be extended to q 2 nax by matching it to a B* dominance model whose normalization is 
fixed with the same sumrule. However, agreement between the various sumrule results is not very good. 

It is interesting to note that HQS, which is so useful for the study of semileptonic B — > £)W decays, where 
it provides a normalization condition at zero recoil and applies to the full, physical, kinematical range, is in 
principle valid for heavy — > light quark decays only around q 2 = q^ax an d provides no normalization. One can 
nevertheless investigate the possibility of using HQS to relate B — » tt£Di decays to D — » tt£Di decays, to obtain 
model- independent information about the former Jl^]-|^]. This approach is taken one step further in Ref. |2j| 
where the structure of the order 1/mi, corrections to the HQS limit description of B — > tt£D^ decays is explored. 

Attempts have also been made to account for the short distance effects which become important at large 
recoils by factorizing B — > n£p£ decays into a hard perturbative part and a soft non-perturbative part, but it is 
unclear whether the b quark is massive enough for this factorization to work f24|| . 

On the lattice one is limited by cutoff effects. Present day lattices have inverse lattice spacings of the order 
of 3 GeV which means that B mesons cannot be simulated directly. One way of preceeding is to perform the 
calculation for heavy quarks with masses around that of the charm and extrapolate the results to rrib using 

2 Here and in the following, £ stands for the e and the fi but not the r. 
3 The Particle Data Group reports 0.002 < \V ub \ < 0.005 [§ 
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heavy-quark scaling laws. One can also try to circumvent this problem altogether by working with discrectizcd 
versions of effective theories such as Non-Relativistic QCD (NRQCD) or Heavy-Quark Effective Theory (HQET) 
in which the mass of the heavy quark is factored out of the dynamics. But even then the momentum of the final 
state hadron in semileptonic B — > it and B —* p decays is large compared to the cutoff in much of phase space 
and, furthermore, these theories are only applicable around q^ax- Thus, either approach will typically yield 
form factors close to q 2 nax and one must perform extrapolations over a large range of q 2 to reach small values 
of q 2 which contribute significantly to the total rates. The problem is even more accute for B — > if* 7 decays 
where one again obtains form factors close to <?^ aa ., q being the photon momentum, while one needs them at 
the on-shell point q 2 = 0. The results of all these extrapolations will depend very strongly on the assumptions 
made about the (^-dependence of the form factors. Assuming nearest pole dominance, ELC and APE |2q] 
have actually calculated the rate for B° — > Tr + £~De decays from the determination of the relevant form factor at 
a single q 2 ~ f 8 — 20GeV 2 . However, their calculations provide no information concerning the <7 2 -dependence 
of the form factor and therefore have large theoretical uncertainties. 

An important step in constraining the g 2 -dependence of form factors in semileptonic heavy — > light quark 
decays was taken by the UKQCD Collaboration (D. R. Burford et al.) in Ref. |^7| who use lattice results for the 
form factors for semileptonic B — > tt£vi and B —>■ K*j decays around q 2 nax , kinematical constraints at q 2 = 
and heavy-quark scaling laws to select amongst different possible functional forms. Though certain forms appear 
to be favored in their analysis of semileptonic B — > ixli>i decays, the situation is less clear for B — > K*-f decays. 
Furthermore, in the case of semileptonic B — » ptvi decays, there are no kinematical constraints at q 2 = for 
the form factors that dominate the rate and their analysis is not possible for these form factors. 

In light of all these limitations and discrepancies, a fully model-independent QCD determination of the q 2 - 
dependence of the form factors for semileptonic heavy — > light decays is very important. To this end we combine, 
in the present paper, QCD dispersion relations with the lattice results of the UKQCD Collaboration (271] and with 
a kinematical constraint and derive model-independent bounds on the form factors relevant for B° — ► ir + t~9t 
decays. The application to semileptonic B — > irtvt decays of the dispersive constraint techniques developed 
in Ref. |29j for semileptonic kaon decays was initiated very recently by C. G. Boyd et al. in Ref. |50|Q. Our 
analysis extends theirs in many ways. The most notable difference is that we use lattice results instead of model 
results to constrain the bounds, thereby obtaining model-independent bounds which can be used to test the 
consistency of experimental results or other theoretical predictions with QCD itself. ^ The difficulty here is 
that one must develop a formalism to take into account the errors on the lattice results. 

We also improve on the results of Ref. |3(J by obtaining bounds for the two form factors required to describe 
B — > Ttlvi decays-not only the dominant one-which enables us to make use of a kinematical constraint to 
constrain the form factors even further. This again requires a generalization of the bounding techniques of 
Ref. ||§. 

Moreover, we decompose the required polarization function according to a more physical helicity basis. This 
again enables us to improve on the bounds of Ref. [[30). And unlike the analysis of Ref. ours does not 
require knowledge of the B*Bir coupling, gB*B-rt- We in fact use our results to put bounds on this coupling. The 
only physical parameters needed in our approach, then, are particle masses and the leptonic B* decay constant, 
/b», also available from the lattice [M. 

The remainder of the paper is organized as follows. In Section || we provide some general background, 
discuss the requirements of HQS and describe the different parametrizations which we later compare with our 

4 Dispersive bound techniques have also recently been applied to semileptonic B — -» D^lv J3T] , [j^ and Aj, — ♦ h c lv p3| decays. 
5 The use lattice results is briefly suggested in the conclusions of Ref. 13. 
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bounds. In Section we describe in some detail the methods we use to obtain dispersive bounds on the form 
factors. In Section |we develop a formalism for taking into account kinematical constraints. As a by-product, 
this formalism enables one to constrain bounds on a form factor with the knowledge that it must lie within 
an interval of values at one or more values of q 2 . In Section ^ we develop a formalism for taking into account 
uncertainties in results used to constrain the bounds and obtain a probability which enables us to dehne bounds 
with unambiguous statistical properties. In Section ^ we combine all of the techniques of the previous sections 
with the lattice results of Ref. ]27j to obtain lattice constrained bounds for f°(q 2 ) and f + (q 2 )- In Section 
we compare these bounds to various parametrizations for the form factors. In Section ||we derive bounds and 
results for the B*Bir coupling. In Section ^ we derive bounds for the total rate for B° — > ir + l~Di decays and 
compare them to predictions of other authors. We also compare some of these other authors' predictions for 
f + {q 2 ) with our bounds on this form factor. In Section [l^ we summarize our main results and discuss how they 
can be improved. 

In Appendix [X], we give explicit expressions for the bounds and derive useful results. These can easily be 
used, in conjunction with the other results in this paper, to obtain model-independent bounds on form factors 
over all of phase space from results obtained by any means in a limited kinematical regime. In Appendix we 
investigate the range of validity of the QCD calculation required to obtain the bounds. We find, for instance, that 
the choice of Q 2 = —q 2 made in Ref. [ j30[ , where q is the momentum flowing through the polarization function 
used to obtain the bounds, leads to uncomfortably large perturbative and non-perturbative corrections. Finally, 
in Appendix ^| we present some of the parameters of the lattice calculation and discuss systematic errors at 
length, including large errors added to cover possible violations of flavor symmetry in the light, active and 
spectator quarks. In phenomenological applications and in situations where systematic errors dominate over 
statistical errors-as they do here-we suggest that they be taken into account earlier than they usually are in 
the analysis of the lattice results, when one controls them better. Though one may loose some information in 
doing so, this approach can also save one from drawing misleading conclusions. Indeed, information about the 
<7 2 -dependence of the form factors f + (q 2 ) and f°(q 2 ) may be lost if some of the systematic errors taken into 
account, before fitting the form factors to different functional forms, are independent of q 2 . | However, since 
the relative proportions of g 2 -independent and g 2 -dependent systematic errors are not known, we believe that 
it is safer to take into account all errors from the beginning. 



2 General Background 

To describe semileptonic B° — > ir + £~i/£ decays one must evaluate the matrix element 

(n + (p')\V»\B°(p)) = (p + p'-q ml q2 ml ) " /+ (q 2 ) + g" ^ ~ ™* f (q 2 ) , (1) 

where V^=u^b and q ~p — p'. Here q 2 runs from rnf epton c^ to Qmax = ( m B — ra^) 2 . f + and f° correspond to 
the exchange of 1~ and + particles, respectively, and must satisfy the kinematical constraint at g 2 =0, 

/+(0)=/°(0). (2) 

For q 2 ~ ^ al , where HQS holds, there are additional constraints. One can easily show that HQS implies the 
following scaling of the form factors with tub at fixed oj — {m 2 B + m 2 — q 2 )/(2mB , m 7V ): 

/V-?LJ = A»s[^ ( 1 + O (*3™)) (3) 
V m B \ V m B J J 

6 I would like to thank J. Nieves for drawing my attention to this point. 
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and 

(4) 



where A°(lu) and A + (oj) are independent of tub- 

Eqs. (0), (|J) and (^) obviously limit the relative functional forms of the two form factors. If we write 

f°(q 2 )=9(q 2 )f + (q 2 ) (5) 

then Eq. @ implies (as long as / + (0) ^ 0) 

5 (0) = 1 (6) 
and the heavy-quark scaling relations of Eqs. (||) and (^) imply the following scaling relation for g: 

Oil -Qmax) = 7+7 x 1+0 • ( 7 ) 

m B A+(uj) \ \ tub J J 

Thus, in choosing parametrizations for f + and f° to compare with our bounds we must ensure that the 
corresponding functions g(q 2 ) satisfy Eq. (^J) and scale like 1/mg as required by Eq. (R). One set of parametriza- 
tions that is consistent with these requirements, that is simple and physically motivated by pole dominance 
ideas, is [£7| 

/V) = //(I - q 2 /™ 2 o) n and f+(q 2 ) = //(l - g 2 /m 2 + )( n+1 ) , (8) 

where n is fixed and / = /°(0) = / + (0). /, m + and m are the parameters to be determined by a simultaneous 
fit to lattice results for f + (q 2 ) and f°{q 2 )- For the scaling relation of Eq. (Q) to be obeyed, m + and m must 
both be equal to ms up to 0(Aqcd/^b) corrections. Moreover, because f°{q 2 ) depends very weakly on q 2 
in the range covered by the lattice results we need only consider the cases n=0 and n—1 which we shall call 
"constant /pole" and "pole/dipole" fits, respectively. 

We will also compare our bounds on / + with the vector dominance form 

which should be an accurate description of the form factor close to zero recoil since the B* pole is so near 
(ms» = 5.32 GeV ^)). Grinstein and Mende ]35| argue, in fact, that this behavior is likely to persist over the 
full kinematical range, having first shown that it does in a combined heavy-quark, chiral and large N c limit. 
That this behavior approximately persists is confirmed by Ball O] who actually calculates the g 2 -dependence 
of / + using three-point function sumrules and also by Belyaev et al. Jig] who obtain this dependence with two- 
point function light-cone sumrules. In both these calculations, however, the mass of the pole is slighly smaller 
than tob*. To find a (^-dependence for f°{q 2 ) which is consistent with the q 2 -dependence of Eq. (JoJ) for f + (q 2 ) 
we note that in the heavy quark limit f~(q 2 ) = —f + {q 2 ), where q 2 f~(q 2 ) = {m 2 ^ ~ m V) (f°(<l 2 ) ~ .f + (l 2 ))- 
This suggests that with f + (q 2 ) parametrized by Eq. (O), f°(q 2 ) may be well described by 



/ + (0) 
1 — g 2 /m| 



/°( g 2 ) = (l-g 2 /m 2 ) i t J"/ 2 , (10) 



where mi and mi are masses on the order of ms ■ Note that Eqs. (|^) and ( p!fj| ) are consistent with the kinematical 
and heavy-quark scaling constraints of Eqs. (j^) and ([?]). In what follows we refer to the combined fit of f + (q 2 ) 
and f°(q 2 ) to the forms of Eqs. (j|) and (^) as the "fixed-pole" fit. 

Finally, we will consider fits to Eqs. ^ and ( |l0| ) for different values of "ms*" to determine the range of 
pole masses allowed by our bounds. 
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3 Dispersive Bounds 



The bounds on f + (q 2 ) and f°(q 2 ) are derived from the two-point function 



II^O) = i J cfx e lq - x {Q\T {V^{x)V^{Q)) |0) 

= -GT^-^On^ + ^n^), (11) 

with — u'y^b and where Ht(V)(q 2 ) corresponds to the propagation of a J p = 1~ (0 + ) particle. In the deep 
euclidean region, i.e. q 2 = —Q 2 with Q 2 3> Aqcd, this two-point function can reliably be evaluated with 
perturbative QCD. The extent to which the condition Q 2 ^> Aqcd can be relaxed is investigated in Appendix 
[b|. The idea, then, is to use dispersion relations to relate this two-point function to the matrix element of 
Eq. (jl|). As long as the indices in EP 1 "^) are treated symmetrically the absorptive part, which is obtained 
by inserting a complete set of hadronic states between the two factors of the vector current, will be a sum of 
positive terms. Thus, the QCD result for this two point function is an upper bound on the contribution from 
the Bit intermediate state. This bound enables us to constrain f + (q 2 ) and f°(q 2 )- 

Let us explore how this works in more detail. The polarization functions Ht,l{(1 2 ) have a branch point 
at q 2 = (rriB + rn^) 2 which marks the Bit threshold. Below this threshold only one pole, corresponding to 
the B* vector meson, is present. Because of its quantum numbers this pole contributes to TlT^q 2 ) but not to 
Hi J (q 2 ). Thus, the spectral functions ImllT ;,^ 2 ), which represent the absorptive parts of ^It,l(i 2 ) an( i which 
are defined by 

-(.9 M V - ft ft) Imn T (? 2 ) + q^q" Imn L (g 2 ) 

= \ J2(^rs (4) (q - pr) (omr)<r|^t| > , (12) 

r 

where the summation extends over all hadron states with the correct quantum numbers, are given by 

ImU L (t) > ~ x ^y/it-t+Xt-t^J^-Oit t+) (13) 



and 

Imllr(t) > 7T 



,3/2 



(.-»»,) 

+ jjx-i- '"-'-^-'-" i/+wr»(t-«+), <i4) 

where t — q 2 and t± = (tub ± tji^) 2 . In Eqs. and ([Ti}) we have limited ourselves to the contribution of the 
B* state with 



2 

777 

(0\V»\B*(r,p)) = -J^e? (15) 

IB* 

and to the contributions of intermediate Btt pairs. These latter contributions appear with an overall factor of 
3/2 because the two states \B°it~) and \B~tt°) contribute in a ratio of 2 to 1 to the spectral functions in the 
limit of exact isospin symmetry. The inclusion of other states would lead to stronger bounds on f + {q 2 ) and 
f°(q 2 ) but would require one to make assumptions concerning the size of their contributions to the spectral 
functions. 

Now, in QCD the polarization functions I1t,l((7 2 ) satisfy the following dispersion relations (q 2 = —Q 2 ): 
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and 

2 



So, combining Eqs. (|^) and @ with Eqs. @ and (|l|), we find 



Xl(Q 2 ) > - I dtk° L (t,Q 2 )\f°(t)\ 2 (18) 



7T 



and 

^ : = -^3 +- / dtk+(t,Q 2 )\f + (t)\ 2 , (19) 
JB* ) [m%, + Q 2 } 71 Jt+ 

where k° L (t) and ki(t) can be obtained easily. 

To obtain the bounds, we proceed along the lines of Ref. ||. Both Eqs. @ and @ give us an upper 
bound on the weighted integral from t = t+ to oo of the squared magnitude of a form factor. Gcncrically, we 
have 

1 r°° 

J{Q 2 )>- dtk(t,Q 2 )\f(t)\ 2 . (20) 
n Jt+ 

To translate this information into a bound on the form factor f(t) for t in the range [0, we map the complex 
t-plane into the unit disc with the following conformal transformation: 



— = J^- ■ (2D 

1-2 Vt+-L V 1 



Then, the integral around the cut in Eq. fl20|) becomes an integral around the unit circle and the physical 
region for semileptonic B — > decays (i.e. t : — > tJ) is mapped onto the segment of the real line 

z : (1 > z G > 0) — > 0, where z Q = z(£=0). Thus, Eq. (|2(i| ) becomes 

J(Q 2 )>/ — |0(z,Q 2 )/(z)| 2 , (22) 

J z|=l ^ 7r * z 

where we have used the fact that k(t, Q 2 ) is a positive definite quantity. Here and in what follows we use the 
freedom we have in defining <f)(z, Q 2 ) to make it real and positive. 
Then we define an inner product on the unit circle via 

(9\h)=! ^g{z)Hz) (23) 

J\z\ = l 27TIZ 

so that the inequality of Eq. ( pp| ) can be written 

J(Q 2 ) > WW) ■ (24) 

The value of the form factor at any point z(t) within the unit circle can be obtained by considering the inner 
product of \(j)f} with the state \g t ) such that 

9t(z) = - \ . (25) 
1 — z\t)z 

If f(t) has no poles below the cut, then Cauchy's theorem yields 

(g t \<Pf) = 0(z(i),Q 2 ) f(z(tj) . (26) 
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Now, because of the positivity of the inner product, we know that 

det 



{4>fW) W\9t) 

(gMf) (gt\gt) 



> 



One can easily show that this condition on the determinant, together with Eq. (|24|), implies 

\f(t)\ 2 < J(Q 2 )-. ' ' 



2 )|- 



(27) 



(28) 



i-K<)r (*(*),< 

which is the bound we are after. The beauty of the method of Ref. [gjl is that it enables one to incorporate 
information about the form factor /(f) to make the bounds more constraining. Suppose, for example, that we 
know the value of /(f) at a discrete set of points t i: i = 1 • • • N. Then, as above, the positivity of the inner 
product guarantees that the determinant of the (N + 2) x (N + 2) matrix 



M(f(t)J) 



( WW) W\9t) W\g tl ) ■■■ W\g tN ) \ 
(gt\(f>f) (gt\g t ) (ftlftO ■■■ (gt\g tN ) 
10/) {gtAgt) (gtAatx) ■■■ (gtAgt N ) 



(29) 



V (gt N \4>f) (gt N \gt) (gt N \g tl ) 



(gt N \gt N ) J 



is positive semi-definite: 



detM(/(f),/)>0, 



(30) 



where / = (/(ii), • •• , /(iw))- This condition on the determinant, together with Eq. (24), leads to an inequality 
for a quadratic polynomial in /(f) which in turn leads to the more restrictive bounds 



Fio(t\f) < /(f) < F up (t\f) , 
where the general expressions for Fi (t\f) and F up (t\f) are given in Appendix [a|. 



(31) 



Now, if f(t) has a pole at t = t p away from the cut, Eq. (|26|) becomes 

(ft = (z(t),Q 2 ) f (z{t)) + K ^f_ Z ^ , (32) 

since a simple pole in f(t) at t = t p translates into a simple pole in f(z) at z = z (t p )^\ The rest of the argument 
leading to Eq. © remains unchanged and we have, instead of Eq. (|3l|), 



Fio(t\f) 



Res ((f> f,z(t p )) 



</»(^(t),Q 2 ) z(t p )-z(t) 



<f(t)< \F up (t\f) 



Res ((/)f, z{t p )) 



0(z(t),Q 2 ) z{t p )-z(t) 



(33) 



where one must not forget the pole contributions to (ftj0/) in the matrix which determines Fi a (t\f) and 
F up (t\f). 

If one does not know the residue of /(t)'s pole at t = t p , one can still obtain bounds on /(f) from the 
knowledge of the position of the pole alone |3(J . All one has to do is perform the replacement 



[z, Q 2 ) -> Mz, Q 2 ) = ^, Q 2 ) " 

1 - Z(tp) z 



(34) 



where t p is assumed to be in the range [f_, f+] so that <p p is positive for z = z(t) with f G [0, f_]. Then, 4> p {z) f{z) 
does not have a pole at z = z(t p ) so that (gi\4> P f) is given by Eq. (pq) with <fi(z) replaced by 4> p {z). Furthermore, 



7 Generalization to the case where f(t) has more than one pole is straighforward. 
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because (z — z(t p ))/(l — z(t p ) z) has magnitude one on the unit circle, (4> p f\(f> p f) is equal to {^>f\4>f) and the 
crucial QCD constraint of Eq. ( |24| ) is left unchanged. Thus, the arguments which lead to the bounds of Eq. (^) 
can be applied here and we find 

F? (t\f) < /(*) < K P (t\f) , (35) 
where Ff (t\f) and F£ p (t\f) are the functions F[ (t\f) and F up (t\f) of Appendix |A] obtained by replacing tj> with 
4> p . Because the bounds obtained with cj> p assume no knowledge of the residue of f(t)'s pole at t = t p , they will 
be looser than bounds which do assume some knowledge of this residue. 

Finally, the bounds of Eq. ( |35| ) can be used to constrain the residue f p of the pole of the form factor f(t) 
at t = t p . All one has to do is multiply the inequalities of Eq. ( |35| ) by 1 — t/t p and take the limit t — > t~ of the 
resulting functions. With the definition 

f[i) = T^Tjr p + fre9{t) ' (36) 



where f re g(t) is finite at t = t p , we find 



with 



fl <f P < P up , (37) 



4 Implementing the Kinematical Constraint 

4.1 The Kinematical Constraint 

In Section || we describe how to obtain independent bounds on the form factors / + (t) and f° (t) given the values 
of these form factors at discrete sets of points. The two form factors, however, are not entirely independent 
since the kinematical constraint of Eq. @ requires that / + (0) = /°(0). In the present section we describe how 
to incorporate this constraint into bounds for f + (t) and f°(t). 

Without this kinematical constraint but given f + (t) at t — t± , ■ • ■ , (/+ = (/+(t+), • ■ • , f + (tZ f+ ))) and 
f°(t) at t = 4,... , t° N0 (/° = (/°(t°), ■ • ■ , f°(t% ))), the methods of Section | yield: | 

F+(t\f+) < f+(t) < F+(t\f+) (39) 

and 

F$,(t\P) < f°(t) < Fl(t\p) . (40) 
Combined with the kinematical constraint, these bounds yield the following for f a = / + (0) = /°(0): 

4>io < fo < 4>u P , (41) 

with 

io = Max j£(t|/°)) (42) 

and 

<f> up = Mm (F+(t\f+),FZ p (t\P)) . (43) 

This is all that we know since the bounds of Eqs. (^) and ( f40| ) carry no indication as to the probability that 
f + (t) or f°(t) take any particular value within them. Thus, implementing the kinematical constraint reduces 
to finding the bounds on f + {t) and f°(t) given /+, f 3 and f a 6 [<pi , (j> U p\- 

8 Here, F stands for either F and F p . 

9 We assume here that Eqs. (^) and ffl) are consistent with the kinematical constraint at t = and will leave to Section |jj the 
discussion of what happens if they are not. 
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4.2 Bounds on a Form Factor Given that It Lies within an Interval at a Discrete 
Set of Points 

To find the bounds on a form factor /(f) given / = (/(fi), • • • , /(fjv)) and Fi a (tN+i\f) < <j>io < x = /(fjv+i) < 
<t>u P < F up (t^ + i\f) (for the case at hand, fjv+i = 0), P| consider the bounds Fi Q (t\f,x) and F up (t\f,x) obtained 
from the constraint detM(/(f), /, x) > using the methods of Section ||. 

According to Appendix [a|, F up (t \f, x) increases monotically with x, for fixed f and /, from x — i^ (fjv+i|/) 
to x = x max where it reaches its maximum value F up (t\f , x max ) — F up (t\f ) and then decreases monotically 
until x = F up (tN + i\f) beyond which it is not defined. Similarly, Fi a (t\f,x) decreases monotically from x — 
Fi (tpf+x\f) to x = x m i n where it reaches its minimum value Fi (t\f,x m i n ) = Fi (t\f) and then increases 
monotically until x = i^ up (fjv+i|/) beyond which it is not defined. 

Therefore, the upper bound on /(f) given / and x £ [<pio, 4>up] may be defined as the maximum of F up (t \f, x) 
w.r.t. x at fixed f. Thus, 

!Fup(t\f 'i 4*up) ^ x raax > 4>up 
F up {t\f) if (44) 

F U p(t\f, (j>lo) if X m ax < 4>lo 

Similarly, the lower bound may be defined as the minimun of Fi a (t\f,x) w.r.t. x: 

{Fio(t\ f ', 4>up) if x m i n ^> (f> up 
Fio{t\f) if (45) 

Fio{t\f,4>io) if x mm < (f>i 

Graphically, these definitions correspond to defining the upper and lower bounds as the upper and lower bound- 
aries of the envelope obtained by considering the set of pairs of bounds (F; (f |/, x), F up (t\f, x)) and allowing x 
to vary freely in the interval [4>i , <j> U p\ ■ 

Furthermore, the behavior in x of the bounds Fi a (t\f,x) and F up (t |/, x) implies that constraining x to lie 
within the interval [4>i ,4>up] can only increase the strength of the bounds, i.e. the kinematical constraint will 
lead to bounds that are at least as constrained as those obtained without it. In practice, we will find that the 
kinematical constraint significantly improves the bounds on / + (f), especially for small f. 

Finally, the conclusions of Appendix |A| can be used to generalize the above results to the situation where 
bounds on a form factor /(f) are sought, given / = (/(fi), • • • , /(fw)) and the constraints (j)\ < /(fi) < (t> l up , 
i = N + 1, • • • , N + M, for any M. 



5 Taking Errors into Account: Bounds and Probabilities 

Up until now, we have assumed that the values / = (/(fi), • • • , /(*/v)) used as input to constrain the bounds 
on form factor /(f) were both exact and consistent with the dispersive constraints. In the present section 
we show how the methods of Section || and Section |] have to be extended to deal with the more realistic 
situation in which the input vector / has errors, as is the case when it is given by the lattice, experiment 
or any other approximate means. Thus, in what follows, we will assume that / is distributed according to 
a normalized probability distribution Pi n (f) [d N f\- When considering the two form factors / + (f) and /°(f), 

10 We choose [4>loj < t > up\ Q [-^iofXiV+l •fuj?(tiV+l 1/5] since, as shown in Appendix |a|, x cannot lie outside 
[FUt N+1 \f),F up {t N+1 \f)]. 
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we will assume that /+ = (/+(£+), • • • , and /° = (f°(tf), ■ ■ ■ , f°(t° N o)) are distributed according to 

V m (f + J )[d N+ f}[d N °f}. 

The problem which arises when / is obtained in an approximation to QCD, from a model or from experiment 
is that the inequality of Eq. (^) may not have a solution. In Appendix |X| we show that a solution exists iff 
detM(/) > 0, where M(/) is the matrix of Eq. ( p9|) with the second row and column deleted, as suggested by 
our notation. That is, a solution exists iff / is itself consistent with the dispersive constraint. The measure for 
/ which ensures such consistency is: 

dfi(f) = [d N f]6(detM(f)) , (46) 

where 9(x] is the standard "theta-funtion" . 

When considering B — + tt£v decays, we must also ensure that the kinematical constraint of Eq. (||) is 
satisfied. The measure for / + and / which incorporates the requirements that / + and f° be consistent with 
the dispersive and kinematical constraints, is: 

dKf + ;P) = [d N+ f + ][d N °f}9(detM + (f + ))e(dctM (p)) 

(F2 p (0\p) - F+(0|/+)) 9 (F+ (0|/+) - F? o {0\p)) , (47) 

where M + (f + ) and M°(P) are the matrix of Eq. ( p9|) relevant for bounds on f + (t) and f°(t), respectively. 

We can use this measure and the probability distribution of the input, Vi n (f + , P) [d N+ f] [d N ° /] , to define 
bounds which have a clear statistical meaning. In doing so, it is important to keep in mind that, for fixed / + 
and p 3 , the bounds given by the methods of Section |^ come in pairs and carry with them no indication as to 
the probability that f + (t) or f°(t) take any particular value within them. 

A probability which takes these two points into account is the probability that a pair of bounds (F^(t\ 
/+,/+(0) e [<t>io,<f>u P ]), F+(t\ /+,/+(0) S [<p lo , 4> UP ])) for f+(t), defined from Eqs. © and © with cp lo and 
<Puv given m Eqs. (H) and (|f), lie within the interval [/+, /+] at a given t. This probability is given by 

V+(t; [ft f+]\f + ; /°) = } / d^(f+; f°) V m (f + , f°) [ftp ~ F+(t\f+, /+(0) G [<f> lo , <f> up ])) 

x (F+(t\f+, f+{0) G [0 io , <j> up ]) - /+) , (48) 

where 

z(f+ ;P) = J d t 4f + ; P) v in (/+ , P) (49) 

is the probability that / + and f° are consistent with the dispersive and kinematical constraints and where 
dfx(f+;P) is given by Eq. @. 

Because only pairs of bounds which lie entirely within [f£ , /+,] are counted, the probability of Eq. ( f48| ) is 
the minimum probability that the form factor f + (t) take a value inside [fj^,fu P ] at t. A similar probability 
can be defined for the bounds on f° (t) and it is these two probabilities that we use in Section || to plot various 
"confidence level" (CL) bounds for /+(*) and f(t). 

6 Bounds on f + (q 2 ) and f°(q 2 ) 

We may now derive bounds for f + (q 2 ) and f°(q 2 )- We proceed in stages to show how the various elements 
that we combine contribute to improving the bounds. Thus, we begin by determining bounds on f + (q 2 ) and 
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f°(q 2 ) without imposing the kinematical constraint of Eq. (0) nor using lattice results. We then impose the 
kincmatical constraint and in a third step only, combine the bounds with the lattice results of the UKQCD 
Collaboration |§7)- 

As described in Section |[ we use the inequality of Eq. ( [Lsl) to obtain bounds on f°{q 2 )- Because the 
Bq(J p = + ) lies above the Bit threshold, its contribution to ImD^i) and f° does not enter the determination 
of the bounds. Thus, without the kinematical constraint nor the lattice results, the bounds are given by Eq. (|2q ) 
with 



where we have used the results of Rcf. |37j an d where (3{t) = \J (t + — t)/(t + — t-). 

The bounds on f + (q 2 ) are obtained from Eq. (|l9|). Combining Eq. ( |l9|) with Eq. ( po|) gives 

1 r°° 

J + (Q 2 )>- dt k+(t,Q 2 )\f+(t)\ 2 , (51) 
* Jt+ 

where 

J + (Q 2 ) = Xt(Q 2 ) - (^V - 2 1 ■ ■ (52) 
V IB- J [m|, + Q 2 } 

We take /b* =33 ^ . in Eq. (|5l|) as given by the lattice in Ref. Q] on the same configurations as those used 
to obtain f+(q 2 ) and f°(q 2 ). 

Now, unlike f°(q 2 ), f + {q 2 ) has a pole at q 2 = m 2 B , < t + . This pole is due to the fact that the virtual 
ly-boson emitted in B° — > ir + £~D£ decays can couple to the Bw vertex through a B*. At leading order in a 
Heavy-Meson Chiral Lagrangian formalism, the B*Bir coupling is of the form |£lj 

C B *Bir=gTr{B a B b Al W } , (53) 

where B a is a field which annihilates pseudoscalar and vector mesons composed of a & quark and a light antiquark 
of flavor a; B a is the dirac conjugate of B a ; and — —d^Tr/f^ + ■ ■ ■ with f n ± = 131 MeV (for details, see 
for example Ref. psfl ). A similar interaction can be written down for mesons composed of b antiquarks and 
light quarks. The contribution of this coupling is straighforward to evaluate. We find that close to g^az ^ 
contributes to f + (q 2 ) a pole of the form of Eq. ( |36| ) with q 2 a i e — m 2 Br and residue 

, IB-B^grh 1 1 

■IP = A> T~ = 9b* B-n^rz , (54) 

V2/7T IB" IjB* 

where fh = JmB^niB and where ^b^+b+tt — — 1 an d 1b*+b°tt+ = 

Since we do not know the value of the coupling g and hence the residue of / + (<? 2 ) at q 2 = m 2 B ,, the bounds 
that do not take into account the kinematical constraint of Eq. (0) nor the results from the lattice are given by 



Eq. (Eq) but with the replacement of Eq. (34) and with 



(!+ 2 ) 2 ( oln s , 1 + z \ 2 for ^ ! + zN 1 



t *{t + -t-)H-z)W\ m + —z) ^- Q) + — Z ) ■ (55) 
0(t) is defined after Eq. @. 

For both f°(q 2 ) and f + (q 2 ), the QCD evaluation of the corresponding subtracted polarization function 
(xl(Q 2 ) and xt{Q 2 ), respectively) is performed in Appendix Because this evaluation appears to break 

llr The error was obtained by combining in quadrature the statistical error with an additional 10% systematic error. 
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Figure 1: Bounds on /°(|(j ,2 |) and f + (q 2 ) versus q 2 with q 2 in GeV 2 . /°(|g 2 |) is plotted to the left of the vertical 
axis with \q 2 \ increasing leftward from to q^ax and f + (q 2 ) to the right with q 2 increasing rightward from to 
q^ax- This way of plotting the bounds is convenient, for it clearly shows whether they satisfy the kinematical 
constraint of Eq. (^J) . The bounds plotted here do not implement this kinematical constraint nor do they take 
into account the lattice results of the UKQCD Collaboration p7| to which we have added systematic errors 
(points) (see Appendix ^|). Errors on these results are smaller than the size of the points. 

down when Q 2 = — 16GeV 2 -the value used in Ref. ]3C|-and because the choice of Q 2 does not change the 
lattice constrained bounds that we derive below significantly, we take Q 2 = where the QCD calculation is 
reliable. 

We plot the resulting, unconstrained bounds for f°(q 2 ) and f + (q 2 ) in Fig. |l|. These bounds are very loose 
and not interesting phenomcnologically. However, since the bounds on f°(q 2 ) are significantly better than those 
on ,f + (q 2 ), it is clear that the kinematical constraint of Eq. (||) will significantly improve the latter. 

To implement this kinematical constraint, we use the results of Section [|. For the case at hand, (f>i — 
F°( g 2 =0|0) and <j) up = F° p (q 2 =0\6), where F°(g 2 =0|0) and F° p (q 2 =0\6) are the bounds for /°(0) obtained 
without any additional input, i.e. those plotted in Fig. [j]. The bounds on f°(q 2 ) are unchanged while those 
on f + (q 2 ) are given by the versions of Eqs. (|44| ) and ( (45[) relevant for f + (q 2 ) with f + = 0. These new bounds 
are plotted in Fig. |^. The kinematical constraint is correctly implemented since the bounds agree at q 2 = 0. 
Furthermore, the bound on f + (q 2 ) is significantly improved, especially at small and intermediate q 2 where phase 
space is largest. Thus, even though f°(q 2 ) does not contribute to the total rate in the limit of vanishing lepton 
mass, its bounds are important because, through the kinematical constraint, they improve the bounds on the 
form factor f + (q 2 ) which determines the rate. 

We now consider the effect of using lattice results to further constrain the bounds. The lattice calculation 
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Figure 2: Bounds on f°(\q 2 \) and f + {q 2 ) versus q 2 with q 2 in GeV 2 , as in Fig. [I]. These bounds implement the 
kinematical constraint of Eq. (||) but do not take into account the lattice results of the UKQCD Collaboration 
to which we have added systematic errors (points) (see Appendix |c|) . Errors on these results are smaller 
than the size of the points. 

of Ref. (27) provides the values of f a (q 2 ) and f + {q 2 ) at three points above q 2 ~ 15 GeV 2 (see Appendix |^). 
Errors are taken into account as described in Section [5] and below. Because of the large systematic errors in 
the lattice results, it is unclear what form the distribution T'(f + ,f ) should take. We have made the rather 
conservative assumption: 1) that the results are gaussian distributed about their central values with a variance 
given by the large errors of Table ^; 2) that these errors are uncorrelated. We have checked that correlations 
have a tendency to reduce the range of the bounds, as one would expect. We have also checked that a step 
function distribution gives more constraining bounds. 

We use the probability of Eq. ( ftjjj ) for / + and the equivalent probability for f° to define "confidence level" 
(CL) bounds for f + {q 2 ) and f°{q 2 )- These p% bounds correspond to pairs of functions (fi^(q 2 ),fu P (q 2 )) an( A 

(f?o(q 2 )JUi 2 )) such that ' for a11 ? 2 > 

V + (q 2 ] {f+(q 2 )J+(q 2 )}\f + ;P)=0Mp. (56) 

To evaluate the various integrals involved, we perform a Monte Carlo in which we generate 4000 independent 
bound samples for f + (q 2 ) and f°(q 2 ). We find that the distributions for the various upper and lower bounds 
at fixed q 2 are more or less bell shaped and nearly symmetric. Since the width of a pair of bounds is typically 
quite small compared to the width of these distributions, the latter are a good guide as to the behavior of the 
probablility of Eq. ( fi"8| ) for /+ and the equivalent probability for f°. Hence, to define our p% bounds, we can 
consider the central p% of these probabilities. We find that the density in bounds increases as p% is varied 
from 95% to about 30%. This is confirmed by Fig. ^ where we plot bounds which have a p% ranging from 
90% for the outermost pair to 30% for the innermost one in decrements of 20%: the space between neighboring 



13 




Figure 3: /°(|<? 2 |) and f + (q 2 ) versus q 2 with q 2 in GeV 2 , as in Fig. |l|. The data points are the lattice results 
of the UKQCD Collaboration |^7j to which we have added numerous systematic errors (see Appendix . The 
pairs of curves are, from the outermost to the innermost, our 90%, 70%, 50% and 30% model-independent QCD 
bounds. These bounds implement the kinematical constraint of Eq. (^) and take the lattice results into account. 
As explained in the text, the fact that they narrow less and less as the percentage is decreased is an indication 
that the most probable behavior for the form factors is within the 30% bounds. 

bounds decreases as p% decreases. Beyond 30%, it is unclear whether this density continues to increase. In 
light of the discussions of Section ||, these p% bounds indicate that there is at least a p% probability that 
f + (q 2 ) lie within [fj^(q 2 ), fi p {q 2 )} at each q 2 and similarly, that there is at least a, p% probability that f°(q 2 ) 
lie within [/^(q 2 ), fu P (l 2 )]- Thus, we conclude that the most probable g 2 -behavior of the form factors is within 
the 30% bounds though, of course, predictions outside these bounds but within the 70% bounds are certainly 
not excluded. 

We have checked that all of our results are stable with respect to the number of Monte Carlo samples. 
We have also constructed lattice-constrained bounds for f + (q 2 ) and f°(q 2 ) that do not take into account the 
kinematical constraint of Eq. (|^) and have found, not surprisingly, that the resulting bounds on f + {q 2 ) are 
significantly worse, especially for small q 2 , while those for f°(q 2 ) are not very different. 

To get some idea of how strongly the bounds depend on the choice of Q 2 , we have also determined the 
bounds for Q 2 = —16 GeV 2 assuming that the MS results of Appendix |b] for xl(Q 2 ) and Xt(Q 2 ) are reliable. 
We find that the interval between the 70% bounds on f°(q 2 ) shrinks by at most 3% for q 2 < 21 GeV 2 but 
shrinks by about 20% around ^ a:I . Similarly, the bounds on f + (q 2 ) shrink by at most 3% for q 2 < 20 GeV 2 
but shrink by about 25% around qm ax - 
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fat type 


£4- /n\ 

/ + (0) 


mg. or m+(GeV) 


mi or m ( GeV) 


X /d.o.f. 


fixed-pole 


0.43 ± 0.04 


m s » 


5.9 ±0.3 


0.4/4 




0.39 ± 0.04 


5.17 


6.2 ± 0.4 


0.2/4 




0.52 ±0.05 


5.9 


5.4 ±0.2 


1.7/4 


cst/pole 


0.48 ±0.05 


5.42 ±0.26 




1.0/4 


pole/dipole 


0.27 ±0.03 


5.79 ±0.19 


6.2 ±0.4 


0.1/3 



Table 1: Fits of the lattice results for f + (q 2 ) and f°(q 2 ) to the various parametrizations described in Section 
|: fixed-pole (Eqs. (|) and ©); cst/pole (Eq. (§) with n = 0); pole/dipole (Eq. (|) with n = 1). All of these 
parametrizations are consistent with the kinematical constraint of Eq. (|J) (i.e. / + (0) = /°(0)) and heavy-quark 
scaling laws. 



7 Comparison of the Bounds with Various Parametrizations 

In Fig. | we plot the 90%, 70% and 30% bounds for f + (q 2 ) and f°{q 2 ) together with fits of the lattice data 
to the "constant/pole", "pole/dipole" and "fixed-pole" parametrizations described in Section^. The results 
of these fits are summarized in Table |l[ The bands which appear in the figure correspond to allowing the fit 
parameters to vary in their 68% CL region. Because of the limited g 2 -range of the lattice results for /°, it 
is impossible to reliably determine the two mass parameters mi and 1712 of the "fixed-pole" parametrization of 
Eq. (0). Therefore, we have allowed mi to vary and have taken mi to be the mass of the nearest pole which 
can contribute to /°, i.e. mi = ms(o+) = 5.46(30) GeV |^6| where the B(0 + ) is the first scalar excitation of the 
B meson. The exact value of mi is not, in fact, very important as the fit parameters of the phenomenologically 
important form factor, f + {q 2 ), depend very little on this value: we have varied mi in the range from ms to 
raj + 1 GeV and have found that the fit to f + (q 2 ) changes by less than one part in a thousand. What is 
important, however, is that with both m\ and mi on the order of mg, the lattice data are fit very nicely and 
the resulting curves lie within the QCD bounds. 

As evidenced by the low values of \ 2 /dof obtained for all three sets of fits, the lattice results alone cannot 
effectively discriminate amongst the different parametrizations, at least when systematic errors are taken into 
account before performing these fits. As mentioned in the Introduction, this early inclusion of systematic errors, 
some of which are independent of q 2 , may be partially responsible for this lack of effectiveness. However, since 
the relative proportions of (^-independent and (? 2 -dependent systematic errors are not known, we have chosen 
to be conservative by including all errors from the very beginning. 

Our lattice constrained bounds do not exclude unambiguously any of the parametrizations either: the best 
fits of all three sets lie within our 70% bounds. However, substantial fractions of the parameter values allowed at 
the 68% CL by the fits are excluded by our 70% and even 90% bounds. Furthermore the best "constant/pole" 
fit lies outside the 30% bounds for a large range of q 2 . In fact, only the best "pole/dipole" fit lies within these 
30% bounds for all q 2 . Since these bounds delimit the region of most probable values for the form factors, 
the "constant/pole" parametrization is least likely to be an accurate description of the form factors while the 
"pole/dipole" form is the most probable parametrization. This is confirmed by the the results of Section || 
where we find that only the best "pole/dipole" fit yields a value of the total rate that is consistent with our 30% 
bounds on that rate: the best "fixed-pole" fit rate only lies within the 50% bounds and the best "constant /pole" 

12 Some of the pole mass values allowed in the 68% CL ellipse for the "constant /pole" fit lead are smaller than g^ aI and therefore 
lead to an f + (q 2 ) which diverges. 
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fit rate, only within the 70% bounds. However, the bounds will have to improve before a firm conclusion as to 
the preferred g 2 -behavior of the form factors can be drawn. 

For completeness, we have also performed a fit of the lattice results to the "fixed-pole" parametrization 
with different values of the pole mass (i.e. tub* in Eq. (j^)). We find that the pole masses allowed range from 
5.17GeV to 5.9GeV. Outside this range, the pole prediction lies either above ("m_B*"< 5.16GeV) or below 
("ms»"> 6.0 GeV) our 70% bounds on f + (q 2 )- We have further allowed the pole mass to vary freely. We find 
/+(0) = 0.35(3), mi = 6.6(6) GeV and "m, B ," = 5.01(14) GeV with a x 2 /dof = 0.004/3. As this fit indicates, 
the lattice data alone may favor a pole mass and a value of / + (0) slightly smaller than those obtained in the 
"fixed-pole" fits of Table [if Thus, f+(q 2 ) could be dominated by the B* pole for q 2 greater than 15 - 20 GeV 2 
but have a q 2 -dependence determined by a smaller pole mass for smaller q 2 , as suggested in Ref. |Q. Such a 
behavior is perfectly compatible with our bounds as are, of course, many other more complicated behaviors. 

We delay making comparisons of our results for / + with those of other authors until Section |^ where we 
also compare predictions for the total rate. In particular, we collect the predictions of other authors for / + (0) 
in Table §. 
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Figure 4: /°(|g ,2 |) and f + (q 2 ) versus q 2 with q 2 in GeV 2 , as in Fig. |l|. The data points are the lattice results 
of the UKQCD Collaboration |27j with additional systematic errors (see Appendix |c]). The pairs of fine solid 
curves are, from the outermost to the innermost, our 90%, 70% and 30% model-independent QCD bounds. 
The thick curves are the results of fits and the bands, the corresponding 68% C.L. regions, for fits to the three 
parametrizations described in Section H (a) "constant /pole" ; (b) "pole/dipole" ; (c) "fixed-pole" with the pole 
mass fixed to tub* ■ 
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8 The Couplings gB*Bn and g 

Having obtained bounds on f + (q 2 ), we can use Eqs. ( j37|) and ( ^8| ) to constrain the B*Bw coupling, gB-B-n- 
Using Eq. (|54| ) and f B * = 33 ^ (see comments after Eq. (|52|)), we find 

- 8.9 ± 1.1 < -6.2 ± 0.8 < -3.4 ± 0.4 < g < 6.7 ± 0.8 < 9.3 ± 1.1 < 11.6 ± 1.4 (57) 

and 

- 724 ± 88 < -505 ± 61 < -275 ± 33 < g B *+B^+ < 544 ± 66 < 751 ± 91 < 942 ± 114 , (58) 

where the outer limits correspond to a CL of 95%, the next interval inward, to a CL of 70% and the smallest 
interval to a CL of 30%. The errors come from the uncertainties on These bounds are weak because we 
do not have lattice results for f + (q 2 ) very near q max - This will be remedied, though, as soon as the lattice 
calculation of Ref. ]27j | is repeated with a sufficient number of light quark masses to enable a numerical chiral 
extrapolation of the form factors (see Appendix |c| for details) and thus provide determinations of the form 
factors closer to q max . 

For completeness, we also give the results for these couplings obtained under the assumption of pure B* 
dominance. As we have seen, B* dominance describes the lattice results well and is consistent with our bounds. 
We find: 

g = 0.35 ± 0.05 (59) 

and 

<7b.+bo»+=28±4. (60) 

Both these results and the bounds include statistical and systematic errors since the lattice results used to 
obtain these results do (see Appendix ^ for details). However, it is important to note that only the bounds 
of Eqs. d57j ) and (|5^) are model-independent and that the errors quoted in Eqs. ( |59| ) and (^) do not include 
the possible deviation of the form factor f + (q 2 ) from pure B*-pole behavior. The results of Eqs. ( |59| ) and ( |60| ) 
must therefore be treated cautiously. 

Nevertheless, the pole result for gB*+B°ir+ compares very favorably with the light-cone sumrule prediction 
of Belyaev et al. ]l6[ ] who find 5s»+so 7r + = 29 ± 3. It is rather high, however, compared to the soft-pion sumrule 
result of Colangelo et al. |39|, gB'+B°T!+ = 15 ± 4, and to the QCD double-moment spectral sumrule result of 
Dosch et al, gB*+B n+ = 14± 4 |4(|. As pointed out in Ref. Q, the sumrule of Ref. can be derived from 
the light-cone sumrule of Ref. [Q. Though the authors of Ref. |L6| and Ref. ^] agree on the sumrule, they 
disagree on its evaluation, the former finding gB*+B°-K+ — 28 ± 6 in the soft-pion limit. For a discussion of this 
disagreement as well as a good summary of recent determinations of gs* b-k-, please see Ref. ]l6[ . 

9 Bounds on the Rate, Comparisons with other Predictions and \ V u b\ 

Our bounds on f + (q 2 ) enable us to constrain the rate dT (B° — > tt + ^vi) jdq 2 . In the limit of vanishing lepton 
mass, which is an excellent approximation for t — e, [i, we have 



where X(q 2 ) = (m 2 B + m 2 — q 2 ) 2 — Am 2 B m 2 and < q 2 < q. 



2 

max ' 
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As was done for f + (q 2 ) and f°(q 2 ) in Section we define the minimum probability that the total rate for 
B° — ► ir^l^&i decays takes a value inside an interval [r; ,r„ p ]. This probability is given by 



V r ([r lo ,T up }\f+;f 



Z(f+;f°) 



<W+; ./*) W + , f) 0[T up - G up ( /'+, /+(0) e <A 



(g ;o (/+,/+ (0) e [<fe,0 up ]) 



(62) 



where 



G„ P (/+,/+(0) € [<fe,0 up ]) 



G 2 F |Kb| 2 
1927r 3 m| Jo 



dq 2 X 3 / 2 (q 2 ) 



Max 



^« + P (9 2 |/ + ,/ + (0) G [<?Wu P ]) , ^o(<7 2 |/ + ,/ + (0) G [fro, <Pup}) 



(63) 



and 



G, o (/+,/+(0)e 



G F \V u b\ 
1927t 3 to| J:i 



d<? 2 A 3 / 2 (<? 2 ) 



0, if F+(<j 2 |/+,/+(0) € [<fe,^ P ]) < < F+(q 2 \f+,f+(0) € [^ ,^ up ]) 



Min 



F+Jq 2 \f + J + (0)& [<f>io,<f>u P ]) , i^(? 2 |/+ / + (0)e[^ o ,«M) 



otherwise, 



(64) 



where F£( - ■ ■) and • •) are the bounds on f + (q 2 ) defined in Section [f| The above expressions for the 

bounds on the rate are involved because the bounds on f + {q 2 ) can change signs with q 2 . 

Now, using the probability of Eq. (|6^), we can define p% CL bounds for the rate. To calculate this probability, 
we use the results of the Monte Carlo for the bounds on the form factors of Section ^. Since, again, the width of 
a pair of bounds is typically quite small compared to the width of the distributions for upper and lower bounds, 
the latter are a good guide as to the behavior of the probablility of Eq. ( |62| ) . We find that these distributions 
are bell shaped, but definitely peaked toward lower values of the rate. Therefore, instead of considering the 
central p% for our defining our bounds, we resort to a different procedure. Starting at the egdes of our collection 
of bounds, we work our way inward toward the most probable pairs of bounds by taking steps of constant size. 
Because the density of bounds increases much faster as we approach the most probable pairs of bounds from 
below than from above, we take the upward step to be much smaller than the downward step so that the 
increase in bound density is more or less equal in either direction. We use 0.2|y u f,| 2 ps _1 upward steps and 
l-OlV^bpps -1 downward steps, because they are a good compomise between size and upward and downward 
balance. We then call p% bound the smallest interval obtained by taking these steps which includes at least p% 
of the bounds. Our results depend little on the exact size of the steps or even on the procedure. For instance, 
taking constant percentage steps of 2% or more on whichever side of the collection of bounds is less dense gives 
very similar results. In all cases, we find that this density increases until about only 30% of the bounds are left. 
Even though the density does continue to increase slightly beyond that point, we prefer to limit ourselves to 
the statement that the most probable values for the rate is somewhere within these 30% bounds. The results 
given by our procedure are well corroborated by the distributions of lower and upper bounds mentioned above. 

We present these results in Table ||. For comparison, we also collect the predictions of various authors for 
the total rate. We specify the g 2 -range for which their calculations are valid and the assumptions made to 
extend the results to all of phase space. We further provide our and their predictions for / + (0). The spread 
of results for the rate is very large, the highest and lowest predictions differing by a factor of about 7. Such a 
spread implies a factor of about 2.6 uncertainty in the determination of \V u b\ and clearly shows the necessity 
for a completely model-independent prediction. 
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V ( R° — > 7T+/~7/^ 


f+(0) 


Details 

L/L UC11J.O 


J. 11J.O VVU1 l\ 


2.4 — > 28 


—0 26 — > 92 


95% CL 




2.8 — > 24 


-0 18^0 85 


90% CL 




3.6 -> 17 


0.00 -> 0.68 


70% CL 




4.4 -> 13 


0.10 -> 0.57 


50% CL 




4.8 -> 10 


0.18 -> 0.49 


30% CL 


QM El (WSB)l 


7.4 ± 1.6 


0.33 ± 0.06 


<? 2 =0; £?* pole 


QM g, (ISGW)] 


2.1 


0.09 


~ cfL,„„\ exponential 


QM 0, (ISGW2)] 


9.6 




~ o 2 ; dinole 


QM 101 




0.26 


o 2 =0 


GM [H 




0.24 - 0.29 


OM+B* Dole 

wivi | i s yj uil 


GM 101 

I II 


9.6-15.2 


0.29 - 0.46 


q 2 =0 ~ 20GeV 2 


SR2 f Toll 


14.5 ± 5.9 


0.4 ± 0.1 




SR2 + 3 ITT1 


4.5-9.0 


27 ± 05 


r/^ — and 

y U emu q. max . 








TnoHifior] R* doIp 


SR3 Ol 

k_/X lu 1U 


3 60 ± 65 


23 ± 02 


« 2 — ~ 5 GeV 2 








+ _B* pole for rate 


SR3 Ol 


5.1 ± 1.1 


26 ± 02 


a 1 — ~ 16 - 20 GeV 2 








-h pole fit w/ 








m™/„=5 25, 10) GeV 


LCSR ifil lfil ITtIi 


8.1 


0.24 - 0.29 


g 2 =0 ~ 15 GeV 2 , then £?* pole 


LAT |27|, (UKQCD)] 1, 


7 ± 1 


0.21 - 0.27 


g 2 ~ 15 ~ 20 GeV 2 








+ "pole/dipole" fit 


LAT H (ELC)] 


9±6 


0.10-0.49 


single point at q 2 ~ 18 GeV 2 








+ pole w/ m po/e =5.29(l) GeV 


LAT g (APE)] 


8±4 


0.23-0.43 


single point at q 2 ~ 20.4 GeV 2 








+ pole w/ m po / e =5.32(l) GeV 



Table 2: Comparison of different theoretical predictions for / + (0) and the total rate for B° — > ir + £~De decays 
(in units of \V u b\ 2 ps" 1 ) . The results in the first row block correspond to our 95%, 90%, 70%, 50% and 30% 
bounds on these quantities as described in the present section and Section [| In column 1, QM stands for 
quark model; SR2 for two-point function QCD spectral sumrules; SR3 for three-point function QCD spectral 
sumrules; LCSR for light-cone sumrules; and LAT stands for lattice. The last column specifies the range of q 2 
for which the calculations are reliable and the assumptions and/or fits made to extend the results to the full 
kincmatical range. 



^ The value for the rate differs from the one quoted in Ref. [g7j . Here, the phase space integral is performed from q 2 = 
to q^n ax = 26.4 GeV 2 assuming flavor symmetry in the light, active and spectator quarks as discussed in Appendix 
Neither the rate nor the value of / + (0) include systematic errors [p8|. 
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Our bounds on / + (0) are not very strong and at the 70% CL they are compatible with all other predictions. 
As discussed in Section |[ the most probable value for / + (0) lies within the 30% CL bounds, i.e. in the interval 
[0.20, 0.48], which is still consistent with most predictions for / + (0), except the ISGW result of Ref. 

Our bounds on the rate, on the other hand, do make some of the other predictions quite unlikely. The 
ISGW model result of Ref. j|] is excluded at the 95% CL; the central value of the sumrule result of Ref. |jl2| , 
at the 70% CL, and the value at the tip of the error bar, at the 50% CL. All other results are compatible with 
our bounds at the 70% CL, including the two lattice results of Ref. ]2j| and Ref. [^6|. They were obtained 
from a single measurement of / + for a q 2 ~ 18 — 20GeV 2 , assuming B* dominance with tub* taken from a 
lattice computation with the same gluon configurations. Though the errors are mainly statistical or due to a 
heavy-quark mass extrapolation of the form factors similar to the one described in Appendix ^|, they were made 
as large as possible to try to accomodate the fact that the g 2 -dependence of the form factor is not determined 
by the calculation. 

Also, for q 2 ^ our bounds on f + {q 2 ) have more to say, some of which could be inferred, in a different 
language, from a x 2 analysis with the lattice data alone. Considering only central values for the predictions 
of other authors, we find, for instance, that the sumrule result Ref. |L4| is excluded by our 70% lower bound 
for over one-third of the kincmatical range from q 2 ~ 12 GeV 2 to q 2 ~ 22 GeV 2 and by the 90% lower bound 
from q 2 ~ 15 GeV 2 to q 2 ~ 20 GeV 2 . The WSB model, or equivalently the lattice result of Ref. @ 0, is 
only excluded by the 70% lower bound for less than a fifth of the kinematical range, between q 2 ~ 15 GeV 2 
and q 2 ~ 21 GeV 2 . The light -cone sumrule result of Ref. Jl6| agrees with our bounds extremely well since 
it very nearly lies within our 30% bounds for f + (q 2 ) for all q 2 . Very good agreement is also found for the 
light-front quark model, non-relativistic hamiltonian (NR) and harmonic oscillator wavefunction (HO) results 
of Ref. H which lie within our 50% bounds for all q 2 . Furthermore, the Godfrey-Isgur (GI) result of Ref. || 
agrees well with our bounds since it very nearly lies within the 70% bounds for all q 2 . However, this agreement 
with the results of Ref. || implies that the B*ir contribution to the S-meson wavefunction, which is neglected 
in Ref. [pj, can only become significant beyond q 2 ~ 20 GeV 2 where our bounds on f + {q 2 ) allow much more 
singular behavior than that given by the results of Ref. || . This observation is corroborated by the light-front 
quark model calculations of Ref. ||], since our bounds on f + (q 2 ) favor the result which relies on the smallest 
value of the coupling g that these authors consider. 

We have also computed the rates obtained from the fits to the lattice data of various parametrizations for 
the form factors performed in Section [?]. The results are summarized in Table ||. As these results show, only 
the "pole/dipole" parametrization gives a central value for the rate which lies within our 30% bounds which, 
together with the observations of Section |7|, indicates that, of the forms we have tried, it appears to be the 
most probable one. Then comes the "fixed-pole" prediction which lies within the 50% bounds and finally the 
"constant-pole" rate which is excluded by the 50% bounds but lies within the 70% bounds. However, once 
again, the bounds will have to improve before a firm conclusion as to the preferred g 2 -behavior of the form 
factors can be drawn. 

Finally, the most probable value for the rate is somewhere in the interval [4.8, 10] | V^t, | 2 ^*s — 1 which corre- 
sponds to the 30% CL bounds. Thus, we can summarize our results for the rate as 

r (B° - w+t-ut) = (4.8 «-> 10) + V -1 (65) 

where the errors are obtained from the 70% CL bounds. Because we have been generous in our estimate of 

13 It is important to note that the actual lattice number from which this result is obtained is compatible, within errors, with those 
of Appendix ^| 

14 These errors are not gaussian. 
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fit type 


1 \L> — ► 7T t 111 unilS OI | y ub \ PS 


cst/pole 


15.2 + 7.5-3.1 


pole/dipole 


8.5 + 3.0-1.9 


fixed-pole 


12.3 + 2.5-2.2 



Table 3: Values for the rate obtained from the fits described in Section [?]. The error bars are 70% errors. 



possible systematic errors, the theoretical error in Eq. ( |65| ) is probably an overestimate and the 60% CL bounds 
or even perhaps the 50% CL bounds may give a better estimate of this theoretical error. For instance, the 
estimate given by the 50% CL bounds is: 

T (B° -> ir+i-ut) = (4.8 ~ 10) + J ° 4 \V ub \ 2 ps- 1 . (66) 

Now, the CLEO Collaboration very recently measured the branching ratio for B° — > tt + £~i , £ decays [Q. 
They found 

B (5° -+ n+i-*) = { (L63 ± °' 46 ± °' 34) X 10 " WSB , (67) 
v ' \ (1.34+ 0.35 + 0.28) x 10~ 4 ISGW 

where experimental efficiencies are determined with the WSB model and the ISGW model ||, respectively. 
Because of the model dependence of their results, a precise statement about \V ub \ is not yet possible. However, 
to illustrate the sort of accuracy our results lead to, we proceed to extract \V U {,\ from these measurements. 
Because our bounds favor very strongly the WSB model predictions for /+ and the rate over those of the ISGW 
model, we consider only the WSB measurement. Using the result of Eq. (|65|), we find 



\V ub \ y/T B o / 1.56 ps = (0.0032 <-> 0.0047) + 0.0007- 0.0007+0.0007 (WSB) , (68) 

where the first set of errors is theoretical (non-gaussian) and the second experimental (statistical and systematic 
combined in quadrature on the average value of \V U b\ given by the 30% CL results). The value of t§o we have 
taken is 1.56(5) ps ]59]. Using the less conservative result of Eq. (pq) for the total rate, we find 



| Kb | \/t S0 /1.56 P s = (0.0032 <-> 0.0047) + 0.0002 - 0.0004 ± 0.0007 (WSB) , (69) 
where the origin of the errors is the same as in Eq. (|68|). Finally, using the 90% CL bounds on the rate, we get 



0.0021 ± 0.0004 < \V ub \ y/T B0 / 1.56 ps < 0.0061 ± 0.0011 (WSB) , (70) 
and using the the 95% CL bounds, 



0.0019 ± 0.0003 <\V ub \ y/T B0 / 1.56 ps < 0.0066 ± 0.0012 (WSB) , (71) 

where the errors shown on the bounds are the experimental errors. 

The determination of \ V ub \ given in Eq. ( |68|) has a theoretical error of 37%, computed by taking the difference 
of the 70% CL bound results over their sum. Though this error is by no means negligible, it is nevertheless 
quite reasonable especially since the result of Eq. (|68| ) is completely model-independent and is obtained from 
lattice data which include a conservative range of systematic errors (please see Appendix ^| for details) . The 
determination of \V ub \ given in Eq. (|69|), obtained under less conservative assumptions (i.e. from the 50% CL 
bounds), has a theoretical error of 27%. 

15 Though this error is not gaussian, it gives a good idea of the sort of accuracy achieved. 
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10 Conclusion 



We have combined lattice results, a kinematical constraint and QCD dispersion relations to derive model- 
independent bounds on the form factors relevant for semileptonic B° — > ir + £~vg decays. Even though the 
contribution of the form factor f° to the rate is negligible, we have found its bounds useful for constraining 
/°(0) = / + (0), thereby making the bounds on f + (q 2 ) stronger. And although our bounds do not unambiguously 
exclude any of the parametrizations for the q 2 -dependence of the form factors that we have tried, [^] they do 
limit the range of possible parameter values and indicate that some of the parametrizations are more probable 
than others. We find, for instance, that the "pole-dipole" parametrization of Section || appears to be more 
likely than the "fixed-pole" parametrization which, in turn, appears more probable that the "constant-pole" 
parametrization. We also find excellent agreement with the light-cone sumrule results of Ref. for f + (q 2 ) 
which almost lie within our 30% bounds for all q 2 . Agreement with our bounds for f + (q 2 ) is also very good for 
the non-relativistic hamiltonian (NR) and the harmonic oscillator wavefunction (HO) results of the light-front 
quark model calculations of Ref. || which both lie within our 50% bounds for all q 2 . Further good agreement 
is found with the Godfrey-Isgur (GI) result of Ref. Q which very nearly lies within our 70% bounds for all q 2 . 
This agreement with the results of Ref. |9| implies that the B*tt contribution to the S-meson wavefunction, 
which is neglected in Ref. ||], can only become significant beyond q 2 ~ 20GeV 2 where our bounds on f + {q 2 ) 
allow much more singular behavior than that given by the results of Ref. . Agreement of other predictions 
with our bounds ranges from slightly less good to much less good. However, the bounds will have to improve 
before a firm conclusion as to the preferred q 2 -behavior of the form factors can be drawn. 

We also use our bounds, as well as the vector dominance results, to determine the B*Btt coupling and 
the coupling g which appears in Heavy-Meson Chiral lagrangians at leading order. While our bounds on 
these couplings are quite weak, our "fixed-pole" fit to the form factors yields g = 0.35 ± 0.05 (Eq. (|5|)) and 
9B*+B°Tr+ = 28 ± 4 (Eq. (|6(])). It is important to note, however, that the "fixed-pole" results rely on the 
assumption that / + is dominated by the B* pole over the full kinematical domain. 

We further derive bounds on the total rate. These bounds disfavor some quark model and sumrule predictions 
with some certainty. At the 70% CL, they also enable the extraction of \ V u b \ from the corresponding experimental 
results with a theoretical error of 37%. Since we have been generous in our estimate of possible systematic errors, 
this 37% is probably an overestimate. It can be reduced to 27% if one is willing to accept the 50% CL bounds as 
a good estimate of the theoretical uncertainties. More generally, this extraction is complementary to the lattice 
determination of \V u b\ from the differential decay rate for B° — > p + l~vi decays suggested in J58|. Furthermore, 
the preliminary determination of the branching ratio for B° — > Tt + £~De decays by the CLEO collaboration Q 
suffers from a rather strong model dependence. Perhaps this model dependence could be reduced by using our 
results. 

As mentioned in the Introduction and described in Appendix |c[ we assume flavor symmetry in the light, 
active and spectator quarks in using the UKQCD Collaboration lattice results of Ref. [^7) for f + (q 2 ) and 
f°(q 2 )- Although there is evidence that such a symmetry is reliable for the values of q 2 we consider, we add 
large systematic errors to the data to cover possible violations. We further add a very large range of errors to 
account for other possible systematic effects. With these systematic errors (and in some cases even without) 
the lattice results of Table [6] are compatible with the very few chirally extrapolated points of Refs. [^5| and J2(| 
within one standard deviation or less. Of course, our whole analysis should be repeated once a complete set of 
reliable chirally extrapolated results are available. 

16 These parametrizations fit the lattice results well and are consistent with heavy-quark scaling laws and the kinematical 
constraint. 
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We add all of these systematic errors to the lattice results before performing our analysis. We believe that 
such a procedure is sensible in phenomenological applications when systematic uncertainties dominate over 
statistical ones. As discussed in the Introduction and in Section |?j this gives a more conservative representation 
of the accuracy of the lattice results. We are thus confident that our final results are correct within these large 
errors. 

We have also investigated the effect of reducing lattice errors. With errors reduced by a factor 4, the 
corresponding lattice constrained bounds would certainly discriminate amongst various parametrizations and 
would enable an extraction of \V u t \ from the experimental branching ratio with a theoretical error of about 20%. 
A more effective means of improving the bounds, however, would be lattice results over a wider range of q 2 . 
This range will increase as simulations on finer lattices become available. Of course, a reliable determination 
of the form factors at or around q 2 = would help tremendously. We refrained, in the present paper, from 
using sumrule or quark model results around q 2 — because the systematic errors of these calculations are so 
different from those of the lattice. We wanted to keep the results model-independent and wished to explore the 
extent to which the lattice alone could constrain B° — > ir + £~i>£ decays. In the future, however, one may want to 
use the methods developed here to combine results for the relevant form factors obtained by different methods 
in different kinematical regimes. 

For the quantities which depend on knowing the (^-dependence of the form factors close to q^a^., such as the 
B* Btt couplings or partially integrated rates above the charm production endpoint, significant improvement in 
accuracy can be achieved by simply repeating the simulation which led to the results of Ref. (27) with a sufficient 
number of light quark masses to enable a numerical chiral extrapolation of the form factors (see Appendix ^ 
for details) and thus provide determinations of the form factors closer to q^^. ■ 

We have also performed an analysis of the QCD corrections-both perturbative and non-perturbative-to the 
relevant polarization functions so as to determine the range of Q 2 in which one may trust the QCD evaluation of 
these functions. We have found that these corrections are under control at Q 2 = but not at Q 2 — — 16 GeV 2 . 

Finally, the methods developed here to take into account the kinematical constraint and errors on the lattice 
results are completely general. They can be used with sumrule or quark model results to extend these results to 
kinematical regimes where these calculations break down. Furthermore, our methods are, in principle, applicable 
to other processes of great physical interest such as the rare decay B — > K*^, where a model-independent guide 
for extrapolating the lattice results for the relevant form factor from q 2 around ^ na . to q 2 = is urgently 
needed, [^] or to the semileptonic decay B — » piv which would enable an independent determination of \V u b\- 
Also, as a by-product of implementing the kinematical constraint, we have derived a formalism which enables 
one to constrain bounds on a form factor with the knowledge that it must lie within an interval of values at 
one or more values of q 2 . More generally, we believe that the combination of lattice results with dispersive 
techniques will lead to many more interesting results in the future. 
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A Explicit Expressions for the Bounding Functions 

Consider the positive semi-definite, (L + 1) x (L + 1) matrix 

/ WW) Wlau) W\gt 2 ) ■■■ W\g tL ) \ 
{9tMf) (gtAgu) (gtAgt 2 ) ■■■ (gtAgt L ) 
M(f)= (gt 2 \(t>f) (gt 2 \gn) (gt 2 \g t2 ) ■■■ (gt 2 \g tL ) 

\ {9t L W) {9t L \9ti) (gt L \gt 2 ) ■■■ (gt L \g tL ) J 



(72) 



with / = (f(ti), • • • , /(tz,)) and all ti distinct. The elements of this matrix can be chosen real for t\,t2, • • • in 
the interval (— oo,t+) in which case the matrix is symmetric. Then, its determinant is 



detM(/) = (0/|0/)detM(/)« 1 ' 1 )> - £ aij (g u \<l>f){g tj \<j>f) 

i,3 = l 



(73) 



with 



dij = (-iy+ 3 det M(/)« 1 - 1) ^ l+1 ' i+1 )> (74) 

where the matrix M(/)^ il Ji')>( i3 '-? 3 )>'"} is the matrix obained by deleting rows «i, *2 ; • •■ and columns ■ ■ ■ 

from M{f). 
Then, 



detM(/) > => J(Q 2 ) det M(/)« 1 ^ 1 )> > ^ ay {guWX^W) 



(75) 



where we have used Eq. (|2~i| ) and the fact that detM(f)^ 1 ' 1 ^ is positive because of the positivity of the inner 
product. 

The positivity of the inner product further implies that: |^] 
1) an > for alH = 1, • • • , L so that Tr a > 0; 



I) det a = (detM(/)« 1 - 1 ») 



L-l 



> 0. 



Therefore, the constraint of Eq. ( f75| ) requires / to lie within the L-dimensional volume in /-space delimited by 
the L — l-dimensional ellipsoid centered at the origin and defined by the equality in Eq. ( f75|) ((j>(z(t), Q 2 ) is 
chosen real and positive for all £ < £+). One can further show that this ellipsoid is circumscribed by the box 
whose sides are given by the hyperplanes which are solutions to the L equations det M(fi) = 0, i = 1, • • • , L. 
What this means is that the maximum and minimum values that any of the f{ti) can take inside this ellipsoid 
are given by the unconstrained bounds of Eq. (^8|) with t = ti. 



Now, suppose that we wish to find the bounds that Eq. ( |75|) imposes on /(ii) given /(£2), • • • , f{th)- Eq. rt7£ 
can be rewritten as 

7^a(. 9tl |#) 2 + 2/3( 5tl |#) (76) 

18 aa ^ and det M(f)^ 1 ' 1 ^ ^ as long as the ti, i = 1, ■ ■ ■ , L are distinct. 
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with 



a = an j 
L 

= ^2 a ii(gu\(t>f) , 

L 

and 7 - J(Q 2 )detM(fy^-J2a l3 (9um(9 tj W)- (77) 

One can show that the discriminant of the quadratic inequality of Eq. ([76]) is 

A = det MC/)^ 1 ' 1 ^ detM(/)« 2 > 2 )> , (78) 
where M(/)« 2 ' 2 )> = M(f(t 2 ), • • • , f(t L )). Since det M(/)« 1 ' 1 )> > 0, Eq. © will have a solution iff (/(t 2 ), • ■ • , 



/(ii)) themselves satisfy the dispersive constraint det M(f(t2), ■ ■ • , /(iz,)) > 0. If this is the case, we find that 

^ (ii|/(i 2 ),- ••,/(*£)) </(*i) <*W(ti|/(* 2 ), ••■,/(*£)) , (79) 

with 



*io/t*(*il/(*a), • • • , /(tx)) = - /H ' , (80) 



since a = an > 0. The bounds of Eq. (]3l]) in Section |3| correspond to those of Eqs. ( [79[ ) and ( pOQ with L = N + 1 
and with the replacements t\ —* t and U — > ti_i, i = 2, • • • , N + 1. 

Now, to impose the kinematical constraint of Eq. (Q) we need to know how the bounds Fi a (t\f,x) and 
F up (t\f, x) behave as a function of x, with / = (f(ti), • • • , /(tw)) fixed and a; = f(tjsr + i). These bounds are 
obtained from the condition det M(f(t), f, x) > 0. The results of this appendix imply that the bounds Fi a (t\f, x) 
and F up (t\f, x) exist iff det M(f, x) > 0. This in turn means that x is constrained to the interval [xi 0l x up ] with 

%lo{up) = Flo(up)(tN+l\f)- 

Furthermore, because the condition det M(f(t), /, x) > implies that (f(t),f,x) must lie on an N + 1- 
dimcnsional ellipsoid, Fi ( up }(t\f,x), viewed as a function of x with t and / fixed, must be the lower (upper) 
segment of an ellipse bisected by the line that goes through the points (xi , Fi a (t\f, xi )=F up (t\f , Xi )) and 
(x up , Fi (t\f, x up )—F up (t\f, x up )). Since the tangents to the ellipse at these points are vertical, F up (t\f, x) must 
increase monotically with x from x — xi to x — x max where it reaches its maximum value and then decrease 
monotically until x — x up beyond which it is not defined. Similarly, Fi (t\f,x) decreases monotically with x 
from x = xi to x — x m i n where it reaches its minimum value and then increases monotically until x — x up 
beyond which it is not defined. Moreover, one can show that F up (t\f, x mux ) = F up (t\f), i.e. the bounds one 
would obtain without any constraints on x — /(f^+i)- Simlarly, Fi (t\f, x m i n ) — Fi (t\f). 



B The Functions Xt,l{Q 2 ) in QCD 

To quantify the O (a s ) corrections to the polarization function, IP 1 ' (q) , we use the results of Ref . |4^, ^3) . 
Because the mass of the b quark, m&, appears already at 1-loop order, the relative size of this leading contribution 
and the O (a s ) correction will depend on the renormalization scheme used for this mass. The calculations of 
Ref. ^] were peformed using an on-shell scheme. To test the reliability of these calculations as Q 2 is 
changed, we perform them in two additional schemes. Thus, the three schemes are: 

19 %min and x m ax may equal x io or x up . 



2G 



1) pole: 



• mb 



pole 



a 8 (fj,) in MS-scheme 



• fi= ^Q 2 + m\ 

2) MS: 

• mj = Wib(fi) with mP h ° le = mb ( 1 + 'Km)) an d 

• a s (^t) in MS-scheme 



In 



3) euclidean (Landau gauge): 

• nib — m b wc (fi) with m b ° le = mi, j 1 

• a s (n) in MS-scheme 



(l+SfM^)) and 00*) = -i 



-hi 



2 | pole 



A* = \/Q 2 + m b is a natural scale for this process and, obviously, <fr = in the pole scheme. 

The non-perturbative corrections, which are due to interactions of the b and u quarks with quark and gluons 
condensates, can be evaluated with the standard diagrammatic techniques of Ref. [Q or with background field 
methods as in Ref. fl43|. Using the diagrammatic techniques, we have checked the results of Ref. fill 45 for 



J. Using the diagrammatic techniques, we have checked the results of Ref. [J43J, 
condensates of dimensions 3 and 4. Combining all of these results and keeping contributions of condensates 
with dimension less or equal to 4, we find (setting m u = 0) 
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xUQ 2 ) 


Q 


OGeV 2 


-16GeV 2 


scheme 


pole 


MS 


Euclid. 


pole 


MS 


Euclid. 


1-loop 


1.3 x 1CT 2 


1.3 x 10 -2 


1.3 x 10~ 2 


2.1 x 10~ 2 


2.1 x 10~ 2 


2.6 x 10~ 2 


0(a a )/l-loop 


16% 


17% 


16% 


35% 


32% 


12% 


total pert. 


1.5 x 10" 2 


1.5 x 10" 2 


1.5 x 10" 2 


2.8 x 10~ 2 


2.8 x 10" 2 


2.9 x 10" 2 


m u (uu) /1-loop 


-1% 


-2% 


-2% 


-11% 


-10% 


-66% 


<a s G 2 )/l-loop 


0.03% 


0.04% 


0.04% 


0.2% 


0.2% 


1.2% 


total 


1.5 x 10" 2 


1.5 x 10" 2 


1.4 x 10~ 2 


2.6 x 10~ 2 


2.6 x 10~ 2 


1.2 x 10~ 2 



Table 4: Perturbative and non-perturbative contributions to the subtracted polarization function xl{Q 2 ) at 
two values of Q 2 . 



, , / x \ 2x , . 2x 
(3 + x) 1 - cc)ln + --In (x) + 5 + 2x+ 
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2^ m 2 + Q 2 x 



(84) 



In Eqs. (82) and (p4|), l(x) is a dilogarithm. 



We have evaluated the expressions of Eqs. ( |8l| ) and ( p3| ) for two values of Q 2 to investigate the range of 
validity of the QCD calculation. One might like to take Q 2 close to the resonance region because the strength 
of the bounds may increase as Q 2 approaches this region as sugg ested by Eqs. ©, @, @ and ©. We 
have chosen Q 2 = 0, which is the value traditionally used when dealing with current correlators involving a 
heavy quark, and Q 2 — — 16 GeV 2 as was suggested in Ref. J3(J. For the masses, we use the results of Ref. 
(rh^mb) = 4.3 GeV) and for the condensates: 



(«u) lGeV = (-0.24 GeV) 3 

n -G 2 ) = 0.02(1) GeV 4 , (85) 



7T 



where the quark condensate is taken from |1J] and the gluon condensate, from |4q| . Our results for the different 
contributions to Xl(Q 2 ) and xt(Q 2 ) are summarized in Table [| and Table ||, respectively. 

At Q 2 = 0, the perturbative and non-perturbative corrections to the one-loop results for xl(Q 2 ) and Xt(Q 2 ) 
are under control, especially when the perturbative MS and euclidean schemes are used. Moreover, there is very 
little scheme dependence in the total perturbative contribution and actual values of xl(0) and xt(0). Thus, the 
QCD calculation appears to be reliable at Q 2 = 0. For Q 2 = —16 GeV 2 , even though the scheme dependence 
of the total perturbative contributions to xl{Q 2 ) and Xt{Q 2 ) is small, the size of the O (a s ) corrections are 
large and it is somewhat doubtful that the perturbative series will converge. Furthermore, the quark condensate 
contributions become quite large and even blow up in the euclidean scheme. So even though the values of Xl(Q 2 ) 
and Xt(Q 2 ) are nearly equal in the pole and MS schemes, one cannot really trust the various expansions at 
Q 2 = — 16 GeV 2 . We therefore restrict ourselves to Q 2 = throughout the paper and use the MS results which 
appear to be the most convergent and which give the loosest constraint on the polarization functions at Q 2 = 0. 



28 





Xt(Q 2 ) 


Q 


OGeV 2 


-16GeV 2 


scheme 


pole 


MS 


Euclid. 


pole 


MS 


Euclid. 


1-loop 


4.2 x 1CT 4 


5.1 x 10~ 4 


5.1 x 10~ 4 


1.0 x 10~ 3 


1.0 x 10~ 3 


2.1 x 10~ 3 


0(a a )/l-loop 


25% 


7% 


6% 


44% 


38% 


-40% 


total pert. 


5.3 x 1CT 4 


5.5 x 10" 4 


5.4 x 10" 4 


1.5 x 10" 3 


1.4 x 10" 3 


1.3 x 10" 3 


m u (uu) /1-loop 


2% 


3% 


3% 


33% 


30% 


352% 


(a s G 2 )/l-loop 


-0.03% 


-0.05% 


-0.05% 


-0.6% 


-0.5% 


-6.2% 


total 


5.3 x 10~ 4 


5.6 x 10~ 4 


5.6 x 10~ 4 


1.8 x 10~ 3 


1.7 x 10~ 3 


8.5 x 10~ 3 



Table 5: Perturbative and non-perturbative contributions to the subtracted polarization function xt(Q 2 ) at 
two values of Q 2 . 



q 2 (GeV 2 ) 


f + (q 2 ) 


f Q (q 2 ) 


14.9 


0.85 ±0.20 


0.46 ±0.10 


17.2 


1.10 ±0.27 


0.49 ±0.10 


20.0 


1.72 ±0.50 


0.56 ±0.12 



Table 6: Lattice results for f + (q 2 ) and f°(q 2 ) at three values of q 2 . Please see text for details on the treatment 
of systematic errors. 



C Lattice Parameters and Discussion of Systematic Errors 

The results for the form factors / + and /° presented in Table [| were obtained by the UKQCD Collaboration from 
60 quenched SU(3) gauge configurations on a 24 3 x 48 lattice at (3 — 6.2 [^7|. The corresponding inverse lattice 
spacing as determined from the string tension is a^ 1 = 2.73(5) GeV p9| . Heavy and light-quark propagators 
were obtained using an O (a)-improved Sheikholeslami-Wohlert action Q and O (a)-improved currents Q. 

Because cutoff effects would be too large, the 6-quark cannot be simulated directly. Therefore, the three-point 
functions used to determine the form factors were calculated at four values of the heavy-quark mass straddling 
that of the charm quark ^] and the corresponding results for / + and f° were extrapolated in heavy-quark mass 
to tub according to Eqs. (||) and (Q) allowing for linear corrections in I/tub and including in the errors the 
variations due to possible 1 / m 2 B corrections J27J . Of the six extrapolated points in Ref . |27j , we have kept only 
those for which the corresponding recoil u> is completely independent of the heavy-quark mass so as to limit the 
introduction of possible systematic effects. |^] 

To keep volume errors under control, one cannot work with arbitrarily light quarks. The strategy, then, 
is to perform the calculation for several values of light-quark mass around the mass of the strange and then 
extrapolate the results to the up and down. However, the limited number of light-quark mass values available 
to the authors of Ref. |^7j made it impossible for them to perform a reliable chiral extrapolation. Therefore, 

20 The corresponding hopping parameters are 0.121, 0.125, 0.129, 0.133. 

21 The value of f + (q 2 ) at q 2 = 20.0 GeV 2 was obtained from a heavy-quark mass extrapolation-identical to the one of Ref. ]27|— of 
the four values of /^(q^ax) corresponding to the four heavy-quark masses already described. These four values of f^iq^ax) were 
obtained, in turn, from a "pole/dipole" fit to the results for f + (q 2 ) and f°(q 2 ) at each individual heavy-quark mass value |28[. 
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the results presented in Table g were obtained with light-quark masses slightly larger than that of the strange 
(flight — 0.14144) and light-flavor symmetry is assumed for both the active and spectator light quarks. 

Dependence of / + and / on light spectator-quark mass should be quite small. This is corroborated by 
the results of J52| for semileptonic D — > K, it decays where the form factors display no statistically significant 
variation as the mass of the spectator cju&rk is changed from slightly above TYistrange 

to zero, while holding 

q 2 fixed. Since it is these very same data-supplemented with results at different values of the charm quark 
mass-which were used in Ref. ]27|] to obtain the lattice results quoted here, we expect spectator-quark flavor 
symmetry to be reasonable. In fact, similar spectator-quark flavor symmetry is observed in the form factors 
which govern B — > K*j decays |53| and in those which govern semileptonic B — > plvt decays where it is 
as small as 5%. We will assume here that uncertainties associated with light-spectator-quark flavor symmetry 
are 10% on both form factors for all q 2 . We do not expect much of a <? 2 -dependence in the corrections to 
spectator-quark flavor symmetry for / + since in a pole dominance scenario, which suits the data quite well, 
the mass of the relevant pole does not depend on the spectator's flavor. For /°, which has a less pronounced 
g 2 -dependence, the situation should be even better. 

For the values of q 2 considered here (see Table ||) , we also expect dependence on the active light-quark mass to 
be reasonable. All recent lattice simulations of semileptonic D decays indicate that the ratio f B ^^{0) / f^,^ K {0) 
is consistent with 1 within errors. For instance, the authors of Ref. j5^] find /c-»jf(0)//D->7r(^) = 1-09 i 0.09. 
Furthermore, the author of the sumrule calculation of Ref. ]54| ] finds that active light-quark flavor symmetry at 
q 2 = is even better for semileptonic B decays: f£^ K (0)/f£^ w (0) = 1.10 ± 0.01 versus fg_ K (0)/fg_^{0) = 
1.01 ± 0.02. And because of the kinematical constraint of Eq. these conclusions also hold for /°(0). To 
get a measure of the size of the uncertainties associated with assuming this flavor symmetry for fixed q 2 away 
from 0, we suppose pole dominance for f + {q 2 )- Such behavior is consistent with the g 2 -dependence of the 
lattice results as well as with our bounds. Then, assuming active-quark flavor symmetry at q 2 — 0, wc find 
that f^iq^/ftattil 2 ) - (1 - q 2 /m 2 B ,)/(l - q 2 /m 2 B ,) ranges from 1.04 to 1.09 for the q 2 in Table | and for 
tub* = 5.42 GeV |!|. Thus, adding 4% to 9% errors to the uncertainty on active-quark flavor symmetry at 
q 2 = should give a reasonable estimate of the possible errors on /+ for q 2 ^ 0. For f°, which has a milder 
<7 2 -dependence than that given by pole dominance, we expect that the dependence on q 2 of uncertainties due to 
the assumption of active-quark flavor symmetry to be mild: they should be of the same order as those on f° (0) . 
Now, to accomodate a possible, slight g 2 -dependence of these uncertainties, as well as possible deviations from 
our pole parametrization for the g 2 -dependence of / + and violations of light, active-quark flavor symmetry at 
q 2 = 0, we will assume that this symmetry holds within 10% at q 2 = 0. Thus, we add in quadrature to the 
uncertainties in / + , errors ranging from 14 to 19% and to those in f°, 10% errors. 

The error associated with ignoring internal quark loops (i.e. quenching) is difficult to quantify. However, 
quenched calculations of form factors for semileptonic D — > K, K* decays (see |52| and references therein) give 
results in agreement with world average experimental values, while quenched calculations of the differential 
decay rate for semileptonic B — > D and B — > D* decays agree very well with experimental data f55| . This 
gives us confidence that errors due to quenching for processes involving heavy-light hadrons are rather small. 
Moreover, the errors due to the uncertainty in the scale discussed in the next paragraph should take some of 
the uncertainties associated with quenching into account since quenching is in large part responsible for this 
scale uncertainty. Nevertheless, to account for quenching errors which are not taken into account by the scale 
uncertainty we add a conservative 10% error to the form factors. 

In obtaining the bounds, we have assumed that the inverse lattice spacing is a -1 = 2.7 GeV, as is given 
by the p-meson mass (a~ 1 (m p ) = 2.7(1) GeV or the string tension (a _1 (cr) = 2.73(1) GeV jH) and in 
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accordance with the choice made in p7| . Other physical quantities lead to slightly different values for the 
inverse lattice spacing, typically in the range 2.5 ~ 2.9 GeV for spectral quantities |5^| . This uncertainty in 
the determination of the scale obviously leads to uncertainties in the lattice results. These errors should be 
relatively small, here, since the form factors are dimensionless quantities. There are, however, two places where 
the scale appears. The first is in expressing the heavy-light pseudoscalar masses in physical units for the heavy- 
quark-mass extrapolations of the form factors. The second is in the determination of q 2 in physical units, but 
only in so far as the scale is used to determine the mass of the final state meson (the "pion"), since we obtain 
q 2 from the dimensionless quantity to through q 2 = m 2 B +m 2 — 2mBm n oj, using the experimental value for rag. 
We have quantified these errors by using the heavy-quark-mass dependence of the form factors given in Ref. |27| 
and the (^-dependences obtained in Section We find that they are less than about 5% for both form factors. 

Finally, there are cutoff corrections of order O (a s amQ) and O (((itoq) 2 ) , where uiq is the mass of the heavy 
quark in the initial meson. For the heavy-quark masses and action used in Ref. pjj ), these lattice artefacts induce 
errors in the form factors which range from 5 to 10% as indicated in the study of Ref. |55| . A systematic study 
of similar effects is also given in |57| . Assuming, for simplicity, that these errors grow linearly with tuq , we have 
estimated the size of the error induced in the value of the form factors at the scale of the b quark. We have 
found them to be about 10%. 

To summarize, the errors on the lattice results for / + and f° given in Table ^| are the quadratic sum of the 
statistical uncertainties given in |27j] ; of a 10% systematic error associated with possible corrections to spectator- 
quark flavor symmetry; of an error ranging from 14% to 19%, depending on q 2 , to account for possible deviations 
from light-active-quark flavor symmetry in /+ and a corresponding 10% error in /°; of a 10% error associated 
with quenching; of a 5% error to take into account uncertainties associated with the determination of the lattice 
scale; and finally, of a 10% error to account for possible cutoff effects. With these systematic errors (and in 
some cases even without) the lattice results of Table ^ are compatible with the very few chirally extrapolated 
points of Refs. [^5) and (26) within one standard deviation or less. Finally, because our errors include both 
statistical and systematic uncertainties, all of the results in the present paper incorporate a full treatment of 
statistical and systematic errors. 



References 

R. Ammar et al. (CLEO Collaboration), CLEO CONF 95-09, EPS0165, July 1995. 
Particle Data Group, Phys. Rev. D50 (1994) 1173. 

N. Isgur, D. Scora, B. Grinstein and M. B. Wise, Phys. Rev. D39 (1989) 799. 
M. Wirbel, B. Stech and M. Bauer, Z. Phys. C29 (1985) 637. 
J.G. Korner and G.A. Schuler, Z. Phys. C38 (1988) 511; 41 (1989) 690(E). 
P.J. O'Donnell, Q.P. Xu and H.K.K. Tung, Phys. Rev. D52 (1995) 3966. 
D. Scora and N. Isgur, Phys. Rev. D52 (1995) 2783. 



C-Y. Cheung, C-W. Hwang and W-M. Zhang, IP-ASTP-16-95 and hcp-ph/9602309 



I.L. Grach, I.M. Narodetskii and S. Simula, INFN-ISS 96/4 and |hcp-ph/9605349 
C.A. Dominguez and N. Paver, Z. Phys. C41 (1988) 217. 



31 



[11] A. A. Ovchinnikov, Phys. Lett. B229 (1989) 127. 
[12] S. Narison, Phys. Lett. B283 (1992) 384 
[13] S. Narison, Phys. Lett. B345 (1995) 166. 
[14] P. Ball, Phys. Rev. D48 (1993) 3190. 

[15] V.M. Belyaev, A. Khodjamirian and R. Riickl, Z. Phys. C60 (1993) 349. 

[16] V.M. Belyaev, V.M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev. D51 (1995) 6177. 

[17] A. Khodjamirian and R. Riickl, talk presented by R. Riickl at "Beauty 95", Oxford, July 1995, to appear 



in the proceedings, MPI-PhT/95-97, LMU 18/95 and [hep-ph/9510294 . 
[18] N. Isgur and M.B. Wise, Phys. Rev. D41 (1990) 151. 
[19] N. Isgur and M.B. Wise, Phys. Rev. D42 (1990) 2388. 
[20] G. Burdman and J.F. Donoghuc, Phys. Lett. B280 (1992) 287. 
[21] M.B. Wise, Phys. Rev. D45 (1992) 2188. 
[22] L. Wolfenstein, Phys. Lett. B291 (1992) 177. 

[23] G. Burdman, Z. Ligeti, M. Neubert and Y. Nir, Phys. Rev. D49 (1994) 2331. 
[24] R. Akhoury, G. Sterman and Y.-P. Yao, Phys. Rev. D50 (1994) 358. 
[25] ELC Collaboration, As. Abada et al, Nucl. Phys. B416 (1994) 675. 

[26] APE Collaboration, C. R. Allton et al., Phys. Lett. B345 (1995) 513; V. Lubicz, private communication. 
[27] UKQCD Collaboration, D. R. Burford et al, Nucl. Phys. B447 (1995) 425. 
[28] D. R. Burford, private communication. 

[29] C. Bourrely, B. Machet and E. de Rafael, Nucl. Phys. B189 (1981) 157. 
[30] C. G. Boyd, B. Grinstein and R. F. Lebed, Phys. Rev. Lett. 74 (1995) 4603. 
[31] C. G. Boyd, B. Grinstein and R. F. Lebed, Nucl. Phys. B461 (1996) 493. 



[32] I. Caprini and M. Neubert, CERN-TH/95-255 and hep-ph/9603414 



[33] C. G. Boyd and R. F. Lebed, UCSD/PTH 95-23 and |hcp-ph/9512363 . 

[34] UKQCD Collaboration, R. M. Baxter et al, Phys. Rev. D49 (1994) 1594. 

[35] B. Grinstein and P. Mende, Nucl. Phys. B425 (1994) 451. 

[36] E. de Rafael and J. Taron, Phys. Rev. D50 (1994) 373. 

[37] S. Okubo and I-Fu Shih, Phys. Rev. D4 (1971) 2020. 

[38] L. Lellouch, E. Sather and L. Randall, Nucl. Phys. B405 (1993) 55. 



32 



[39] R. Colangclo, F. De Fazio and G. Nardulli, Phys. Lett. B334 (1994) 175. 
[40] H. G. Dosch and S. Narison, Phys. Lett. B368 (1996) 163. 

[41] L. Lcllouch in the proceedings of the XXXIV th Cracow School of Thcorectical Physics, Acta Physica 
Polonica B25 (1994) 1679. 

[42] S. C. Generalis, Open Univserity preprint OUT-4102-13. 

[43] L.J. Reinders, S. Yazaki and H.R. Rubinstein, Phys. Lett. B103 (1981) 63. 

[44] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385. 

[45] S. C. Generalis, J. Phys. G16 (1990) 367. 

[46] S. Narison, Acta Physica Polonica B26 (1995) 687. 

[47] C. Becchi, S. Narison, E. de Rafael and F. J. Yndurain, Z. Phys. C8 (1981) 335. 
[48] E. Braaten, S. Narison and A. Pich, Nucl. Phys. B373 (1992) 581. 
[49] UKQCD Collaboration, C. Allton et al, Nucl. Phys. B407 (1993) 331. 
[50] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572. 

[51] G. Heatlie, G. Martinelli, C. Pittori, G. C. Rossi and C. T. Sachrajda, Nucl. Phys. B352 (1991) 266. 
[52] UKQCD Collaboration, K. C. Bowler et al, Phys. Rev. D51 (1995) 4905. 
[53] UKQCD Collaboration, K. C. Bowler et al, Phys. Rev. D51 (1995) 4955. 
[54] S. Narsion, Phys. Lett. B337 (1994) 163. 

[55] UKQCD Collaboration, K. C. Bowler et al., Phys. Rev. D52 (1995) 5067; L. Lellouch (for the UKQCD 
Collaboration) in the proceedings of the 30 th Recontres de Moriond: QCD and High Energy Hadronic 



Interactions, Meribel les Allues, France, 19-25 March 1995, Marseille CPT-95/P.3196 and frep-ph/9505423 
[56] UKQCD Collaboration, C. Allton et al, Phys. Rev. D49 (1994) 474. 

[57] M. Crisafulli, V. Lubicz, G. Martinelli and A. Vladikas, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 400. 
[58] UKQCD Collaboration, J. M. Flynn et al., Nucl. Phys. B461 (1996) 327. 

[59] I.J. Kroll, FERMILAB-CONF-96-032, to appear in Proceedings of the XVII International Symposium on 
Lepton-Photon Interactions, Beijing, China, 10-15 August 1995. 



33 



